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ABSTRACT 

We propose a model of allocating galaxies in cosmological N-body simulations. We identify each 
subhalo with a galaxy, and assign luminosity and morphological type assuming that the galaxy lumi- 
nosity is a monotonic function of its host subhalo mass. The morphology assignment is made by using 
two simple relations between subhalo mass and galaxy luminosity of different types. One is using a 
constant ratio in luminosity of early (E/SO) and late (S/Irr) type galaxies at a fixed subhalo mass. 
And the other assumes that galaxies of different morphological types but having an equal luminosity 
have a constant ratio in their subhalo masses. We made a series of comparisons of the properties 
of these simulated galaxies with those of the SDSS galaxies. The resulting simulated galaxy sample 
is found to successfully reproduce the observed local number density distribution except for in high 
density regions. The luminosity function is studied as a function of local density. It was found that 
the observed luminosity functions in different local density environments are overall well-reproduced 
by the simulated galaxies. Discrepancy is found at the bright end of the luminosity function of early 
types in the underdense regions and at the faint end of both morphological types in very high density 
regions. A significant fraction of the observed early type galaxies in voids seems to have undergone 
a relatively recent star formation and became brighter. The lack of faint simulated galaxies in dense 
regions may be due to the strong tidal force of the central halo which destroys less massive satellite 
subhalos around in the simulation. The mass-to-light ratio is found to depend on the local density in 
the way similar to that observed in the SDSS sample. We have found an impressive agreement between 
out simulated galaxies and the SDSS galaxies in the dependence of the central velocity dispersion on 
the local density and luminosity. 

Subject headings: Cosmology: simulation: halo-galaxy: luminosity function; Numerical 
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1. INTRODUCTION 

The current galaxy formation paradigm can be char- 
acterized by "hierarchical clustering" . This means that 
massive dark matter halos form by merging less massive 
halos and/or by accreting ambient matter, and also that 
a dark matter halo governs the evolution of the galaxy 
residing inside. Most galaxies are believed to be hosted 
by the dark matter halos becuase the halos can provide a 
deep potential well for baryonic matter to condense and 
cool down. This leads to triggering of the star formation; 
gas sufficiently accumulated in the halo potential cen- 
ter begins to experience many hydrodynamic processes, 
such as radiative cooling, star-formation, supernova ex- 
plosion, and chemical enrichments. All of these processes 
play an important role in making the visible galaxies in 
the end. Because galaxies as a building block of the 
large scale structures consequently follows the evolution 
of their host halos over the cosmic history, understanding 
gravitational evolution of dark halos is very important for 
the study of galaxy formation and evolution. 

Over the past few decades, cosmological simulations 
have been proved useful for the study of structure for- 
mation. Simulations have boosted up many investi- 
gations of the nonlinear structure evolutions and their 
results have been extensively compared to the obser- 
vations in various aspects of interest. Many studies 
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have reported successful recovery of observational fea- 
tures of the galaxy distribution like the two-point cor- 
relati ons of galaxies (IConrov et al.l [20061 : iKraytsov et al 
2004 [B crlind & Weinberg 2002), topologv (|Park et al 



2005al a lGott et al.ll2006l) . and the environmen tal dcpcn 



dence of spin distributions (jCervantes-sodi et~a l. 2007). 

Detailed modeling of galaxy evolution and understand- 
ing of the galaxy properties have been possible with 
the recent advent of the huge redshift surveys such as 
Sloan Digital Sky Survey^ (SDSS) and 2dFGRS6. These 
larger surveys provide a rich information on the for- 
mation and evolution of galaxy. To meet a require- 
ment to establish a relation between simulated structures 
and observed galaxies many techniques have been intro- 
duced. The semi-analytic model(SAM) of galax y for- 
matio n (jCole et al.lll994l : iKauffmann et al. 1997; Baughl 
l2006( l is based on the hierarchical clustering of halos 
whose merging history trees are built by generating a 
set of Gaussian random numbers. The mass growth his- 
tory of halos on various mass scales can be generated and 
traced by this method. Numerous SAM parameters are 
implemented in the merging history to reflect the hydro- 
dynamical processes of heating, cooling, star formation, 
and aging the stellar populations. Their parameter val- 
ues are fine tunned for the best description of collective 
properties of observed galaxies. But as the number of 
observational contraints increases, the number of param- 
eters increases and the SAM becomes much complicated. 
Hydrosimulations employ direct formulations of the 

* http://www.sdss.org 

^ http://www.aao.gov.au/2df/ 
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hydrodynamic processes (jWeinberg et al.1 [2006? ). Two 
types of hyd ro simulations are now widely adopted; the 
Lagrangian (|Monaghanl |1992|: |H ernquist fc Katz) 119891) 
and E ulerian methods ( Tasker fc Bryan 2006; Hart en. A1 
19971 and for comparati ve study, see iHcitman n et al.l 
2005t iThacker et al.l[2000| ) . They set gas particles or gas 



grids in the system of interest. Then, gas dynamics are 
taken into account by calculating the hydrodynamic in- 
teractions between neighbor particles or by solving dif- 
ferential hydro equations between adjacent grids. This 
method has several weak points; it suffers from the lack of 
resolutions in space and mass as the N-body simulations. 
And in most cases, the star formation and supernova ex- 
plosions are far beyond the simulation resolution. Some 
processes such as radiative transfer, star formation, su- 
pernovae feedback, and initial stellar mass function, are 
not thoroughly understood and are difficult to parame- 
terize in a well established way. 

One of the offsprings of t he SAM is the Halo Oc- 
cupation Distribution (HOD; Z heng et all 120051 ; iSeliak I 
[2000; Berlind fcJVeinberg 2002) model. It relates the 
Friend-of-Friend (FoF) halo mass to the number of con- 
tained subhalos (or galaxies) using a conditional proba- 
bility, P{N\M), where M is the FoF halo mass and N 
is the number of subhalos inside the FoF halo. This 
probabi hty function is obtained from n umerical simu- 
lations (iKravtsov et al.l l2004 iBerlind fc W einberg 2002; 
Jing et al.l Il998l) or f rom observations (.Zehavi et al.l 



20041; lAbazajian et al.l l2005t IZheng et all |2007D fitting 



the model to the observed two-point correlation func- 
tions. 

Another variant of the SAM is a one that 
adopts the one-to-one monotonic correspondence be- 
tween the galaxy luminosity and subhalo mass 
jMarinoni fc Hudsod Hool; I Vale fc Ostrikerl [200l . [20061 ; 
[Shankar et al. 200f|). This method is often called the 
"correspondence model" . It is simpler than other meth- 
ods because it requires only two prerequisites; the lumi- 
nosity function of galaxies and mass function of subhalos. 
Its key assumption that a more massive subhalo hosts a 
brighter galaxy, is consistent with the hierarchical clus- 
tering picture. As the halo mass grows, its luminosity 
is expected to grow in general because merging of halos 
may be followed by the merging of galaxies. By match- 
ing these two functions, the subhalo mass is mapped to 
galaxy luminosity. 

However, sometimes a galaxy may survive the merger 
event because baryonic component of a galaxy is usually 
more concentrated and is more tightly boun d than its 
dark counterpart. These "orphan" galaxies (iGao et ahl 
|2004[ ) which have no separate corresponding halos could 
exist in the cluster regions. Also halos of mass below a 
certain characteristic scale are unable to provide baryonic 
matter with gravitational attraction strong enough to re- 
sist against supernova explosions which can tear up the 
small-mass system. Even though such cases are many, 
it is sound to expect that the galaxy census is closely 
related to the population of subhalos because most of 
the observed galaxies are field galaxies who have their 
own dark halos. It is worth investigating the hypothesis 
of the subhalo-to-galaxy correspondence model of galaxy 
formation. 

In this paper, we apply the subhalo-galaxy correspon- 
dence model to a large N-body simulation and compare 



the statistical properties of simulated galaxies with those 
of galaxies observed by the SDSS. We organize this paper 
as follows. In section 2 we describe our simulation and 
the subhalo finding method. In section 3, we implement 
the subhalo-galaxy correspondence model and show how 
to assign morphological types to mock galaxies. In sec- 
tion 4, the local density is introduced to quantify local 
environments and the local density distributions of mock 
and SDSS galaxies are investigated. We also compare the 
luminosity distribution of simulated galaxies with those 
of the SDSS galaxies in section 5. The environmental 
dependence of central velocity distributions are studied 
for the mock and SDSS galaxies in section 6. And dis- 
cussions and conclusions follow in section 7. 

2. SIMULATION AND HALO FINDING 

We have carried out a cosmological N-body simula- 
tion of the universe with the WMAP 3-year cosmolog- 
ical parameters. We ran 1024'^ cold dark matter par- 
ticles using an impro ved version of the GOTPM code 
(jDubinski et al.|[2004D . The code adopts a dynamical 
domain decomposition for the Particle-Mesh part with a 
variable width of the z-directional domain slabs. It also 
uses more compact and efficient oct-sibling tree walks re- 
ducing the computational cost for the short-range force 
update which consumes about 90 % of total run time. 
As a results, this new version outperforms the previous 
version by about factor three in speed. The simulation 
was run on a beowulf-type system installed at Korea In- 
stitute for Advanced Study. The linux cluster consists of 
256 AMD cores and 1 tera-byte main memory. 

The simulation and cosmological parameters adopted 
in this study are lists in Table [TJ The number of time 
steps in the simulations are empirically pre-determined 
to satisfy the requirement that the maximum displace- 
ment of particles in a step should be less than the force 
resolutions which is 0.1 times the mean interparticle sep- 
aration (dmcan)- The Starting epochs of simulations are 
chosen to constrain that any particle should not over- 
shoot neighbors when the Zel'dovich displacement is 
made. 

Subhalos are i dentified by the PSB method 
(|Kim fc ParkI I2006D . The method apphes the FoF 
algorithm to identify dark matter particle groups adopt- 
ing the standard linking length, /link = 0.2(iinean- Then it 
divides each FoF halo into subhalos. This two-sta ge halo 
finding is in comrnon wi th other methods ( Springcl et al.l 
120011: IShaw et al.ll2006f ). In the second stage we build 
a particle density field by a coordinate-free method as 
usually adopted in the smoothed particle hydrodynamics 
((Monaghan 1992). This adaptive kernel implementation 
is intended to resolve tight clusterings of particles. Then, 
we search for 26 nearest neighbors at each particle. To 
construct isodensity contour we move along the neighbor 
positions in a similar way used in the origingal PSB 
method (for details see Ki m fc Park! 2006). Self-bound 
and tidally stable subhalos are identified by measuring 
the tidal radius of the subhalos and total energies of 
their member particles; member candidates are selected 
if their distances to the center of a subhalo are less 
than tidal radius of the subhalo. Among these particles, 
gravitationally unbound particles are discarded. The 
resulting subhalos are used to allocate galaxies. We 
call the most massive subhalo a central halo and other 
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subhalos satellite halos in the FoF group. Also, the 
group of particles found by the FoF method are named 
the FoF group or FoF halo. 

3. THE SUBHALO-GALAXY CORRESPONDECE MODEL 

We use the monotonia one-to-one correspondence 
model between galaxies and subhalos; there is one and 
only one galaxy in each subhalo and a more massive sub- 
halo hosts a more luminous galaxy. We apply this model 
to assign galaxies within our simulation box. 

The one-to-one correspondence model is formalized in 
the following way. We first measure the mass function 
$(Mft) of FSB subhalos. We then take the observed 
luminosity functions of the early (E/SO) and late (S/Irr) 
type galaxies. The relation between the halo mass and 
the corresponding absolute magnitude limits is given by 
the following integral equation. 

^{M')dM'= / (t>E{M')dM' 



ML{Mh) 



(l)L{M')dM', (1) 



where $ is the subhalo mass function and c/je and 0l are 
the luminosity functions of early and late type galaxies. 
In this equation we separate the full luminosity function 
of subhalos into those of early and late types to take into 
account the difference in mass of the halos associated 
with galaxies of different morphological type. 

Two models are proposed to utilize the simple mass 
and luminosity scaling ratios between different morpho- 
logical types at the fixed luminosity and at the fixed 
mass, respectively. The first model uses 

LLiMh) = Le{^Mh), (2) 

where Alh is the subhalo mass, and and Le are the 
luminosities of late and early types, respectively. The 
equation means that at an equal luminosity the early 
type galaxy is k times more massive than the late type 
galaxy. This model is based on the findings of Park et 
al. (2007a) who have assured that the early type galaxies 
brighter than Mr = —19.5 have about \/2 times higher 
central velocity dispersion and pairwise peculiar velocity 
difference compared to the late types of the same bright- 
ness. This implies k « 2 if {Mh/L)E = K{Mh/L)L. An- 
other convencin g observational evidenc e for k = 2 can be 
also found in Mand elbaum et all (2006) who obtained the 
same result for halos of mass greater than 10^^/i~^Mq 
from the analysis of the galaxy-galaxy weak lensing of 
the SDSS sample. Equation ^ then makes equation ^ 
a one-to-one relation between halo mass and galaxy lumi- 
nosity. From now on, we name it k model. The second 
model uses a constant factor for the luminosity rather 
than the mass. This relation can be formulated as 

LL{Mh) = PLsiMh), (3) 

which says that late type galaxies are brighter than early 
type galaxies of the same mass by a factor of (3. Now we 
call it the /3 model. 

For the galaxy luminosity distribution, we use the 
Schechter function 

(j)iM) - (0.4 In 10) 0*iO-° -*(-^--^*)("+i) 

exp \^io-OAiM-M*) 



where we adopt the type-specific Schechter function pa- 
rameters Me - Slogio^ = -20.23, (1)^ = 7.11 x 
10~'^ft.'^Mpc~"^, and as = —0.53 for early type galax- 
ies and Ml - 51ogio/i = -20.12, (j>l = 12.27 x 
IQ-^/i^Mpc"^, and ul = -0.90 for late types in the 
r band, which are given in Table 2 of Choi et al. (2007) 
for the SDSS galaxies brighter than Mr = —19.0 (Here- 
after, we will drop the term blogigh in the absolute 
magnitude and the magnitudes are all r-band magni- 
tude.) These parameter valu e s are quite different from 
those given bv iBlanton et al.) ()2001l ) and Blanton et al. 
(2003) mostly because we are using the type-specific 
functions. Reducing the effects of the internal extinc- 
tion by Choi et al. also makes a significant difference 
in the a parameter compared to Blanton et al.'s results. 
Difference also comes from the data size (Dr4plus versus 
~DR1), adopted cosmology {Q„i = 0.27, Oa = 0.73 ver- 
sus flrn = 0.3, r^A = 0.7), and sample definition (volume- 
limited versus apparent magnitude-limited). 

The /3 model can be easily solved by using the rela- 
tion, MlIMh) = ME{Mh) - 2.51ogio/3 in equation 
but the K model is a bit tricky to solve and needs a few 
assumptions. The k model can be solved by a chain of 
equations whose n-th term is. 



^{M)dM-^ 



(j)LiM)dM 



(f>EiM)dM, 



(5) 



where Mn-2 < Mn-i < Mn and 

Mh{ML,n-2)=Mh{ME,n-l) = m 

Mh{ML,n-i) = Mh{ME,n) = tti/k, 



(4) 



(6) 
(7) 

for K > 1. Here Mh{ML) and Mh{ME) denote the halo 
masses of late type and early type galaxies of absolute 
magnitude M, respectively. If Mh, Mn-2, and Mn-i 
are given, then Mn can be derived. Figure [1] depicts the 
relations among the three magnitudes described above. 
The solid and dashed curves show the luminosity func- 
tions of early and late types, respectively. And the dot- 
ted line connects the same mass of the two morphological 
types. The shaded areas denoted by An and B„ are the 
integrated number densities of galaxies of magnitude be- 
tween Mn-2 and Mn~i for the late type and between 
Mn-i and Mn for the early type, respectively. There- 
fore, the sum of the two shaded areas should be equal to 
the integrals of the subhalo mass function (the left-hand 
side of Eq. [5]). By making a stride to the next mass 
scale which is smaller by 1/k times, we can get a chain 
of equations and derive the scaling relation between the 
mass and the luminosity. To solve this chain of equations 
we set initial conditions under the plausible assumption 
that the late-type contributions to the number density 
at the high-mass (or bright) end are negligible compared 
to those of early types (or ~ Ai,A2 ^ Bi,B2). This 
setting is quite fair from the fact that early type galaxies 
dominate the bright end population of galaxies. Then, 
we are able to solve the series of equation ([S]) from the 
initial conditions where only the early-type contribution 
to number density is dominant. 

Now we investigate the early type fractions as a func- 
tion of the host subhalo mass. The probability of a 
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TABLE 1 
Simulation parameters 



name 


Np 




Nstep Zi 


h 




nt 




b nip'' e<i 


H2 


1 10243 


10243 256 


3800 95 


0.732 


0.238 


0.042 


0.762 0.958 


1.314 1.0 X 10** 25/i-ikpc 



in h~^Mpc 
^ spectral power index 
particle mass in 
force resolution 




Fig. 1. — Illustration to derive the M^-AIr relation for the k. 
The solid and dashed curves are the luminosity functions of early 
and late type galaxies, respectively. And the dotted line connects 
the same mass of the two types. The shaded areas tagged by An 
and B„ are the numbers of galaxies in the magnitude intervals 
between M„—2 and Mn—i for the late type and between Mn-i 
and M„ for the early type, respectively. 



galaxy of mass to be early type is given by 

(t)EiME\Mh) 



fE{Mu) 



(8) 



where the denominator is the luminosity function of all 
galaxies of mass Mh and the numerator is the luminosity 
function of early types of mass Mh- Figure [5] shows the 
results for the k and /3 models. Using Je we are able to 
randomly assign morphological type to each mock galaxy 
in accordance with its halo mass. For a larger value of /3, 
the early type galaxy tends to dominate the population 
down to lower mass scales (Mc) and below Mc Je drops 
more rapidly. But for the n model, those changes in Je 
are not so much steep. Also, in this plot, we note that a 
more dramatic change of the fE in the [3 model than the 
K, model. 

Figure [3] shows the relation between the subhalo mass 
Mh and absolute magnitude of early [solid lines) and 
late [dashed lines) type galaxies derived from equation 
([1]) and the luminosity functions of the SDSS galaxies in 
the n [upper) and [3 [lower panels) models. Only shown 
are the cases oi (3 — 2 (thick lines) and 1.5 (thin lines) and 
K = 2 (thick lines) and 3 (thin lines). For a comparison, 
the characteristic minimum mass of the central subhalos 
at each absolute magnitude limit estimated by Zheng 
et al. (2007) based on the HOD model is given as filled 
circles. They obtained the magnitude-to-mass relation of 
the central subhalo but ignored the morphological types 
of galaxies. The HOD model is quite consistent with our 
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Fig. 2. — Early type fractions as a function of host subhalo 
mass. Each curve shows the fractional distributions of early type 
galaxies for /3 = 1, 1.5, 2, and 3 {lower panel) and for k = 1, 1.5, 
2, and 3 {upper panel) models. The model parameter values are 
attached to the curves. 

Mh-Ai relation for the early type galaxies with (3 = 2, 
and K = 2 and 3. 

In the lower panel of Figure [4] shows the mass versus 
luminosity relation (solid line) for the earl y type galaxies. 
Here we adopt M.q = 4.64 — Slog^Q h (Bl anton fc Roweil 
l2006f ) to transform the magnitude to luminosity in the r- 
band filter. It can be noted that galaxy luminosity drops 
steeply at the low mass end and rises in a power law at 
the high mass end. Taking into account these features, 
we propose a function 



L[M) 



M* 



ml 



M* 

IT 



7-1 



-[^) 



(9) 



as a fitting function for the halo mass versus early-type 
galaxy luminosity relation. We apply the x-square fit- 
ting to the relation and obtain the best-fit parameter val- 
ues of = 38.3/iMq/Lq, M* = 2.04 x W^^H-^Mq, 
and 7 = 0.644 for the early type in the (3 = 2 model. 
In the K = 2 model, ^rni = 39.8(20.3)/iMq/Lo, M* = 
4.83(1.52) X lO"/i-iM0, and 7 = 0.667(0.719) for early 
(late) type. Only the fitting result for the early type 
in the (3 = 2 model is shown by a short-dashed line. 
The top panel shows the mass-to-light ratios of the early 
(solid curve) types in the (3 = 2 model. The long-dashed 
and dotted curves show those ratios of early and late 
type galaxies in the k = 2 model. In the (3 = 2 model, 
there is a upturn of the mass-to-light ratio of early types 
around M„ = 3 x 10^^/t~^M(Ti (for those measured in 
other bands, see Ivan den Bosch et "all 120031 : lYang et al] 
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Zheng et al (2007) 



Early(/?=2) 

Edrly(^ = 1.5) 

Lale(fl = 2) 

■ Late(ff=l,5) 



10* ^ 



1012 10'= lO" lO's 

Fig. 3. — Scaling relation between halo mass and galaxy mag- 
nitude in the ft {upper panel) and /3 [lower panel) models. In the 
lower panel, solid lines shows the dependence of the modeled early 
type magnitude on the host subhalo mass for (3 = 2 (thick) and 
for /3 = 1.5 (thin). And in the upper panel, we show the mass- 
to-magnitude relations for k = 2 (thick) and k = 3 (thin) models. 
Dashed curves show the relations of late type galaxies. Filled cir- 
cles are the best-fitting results of the SDSS observations to the 
HOD model in iZheng e t al. (2007). And the errorbar is the width 
of the cut-off profile of the step-like function in the HOD model 
(for details, please see their paper). 

[200aiEke et"all[200l . which means the star formation 
in galaxies is strongest at the mass scale. A fitting for- 
mula for the mass-to-light ratio is derived from equation 

^Ve^^)- (10) 

Using this equation, the mass scale corresponding to the 
minimum mass-to-light ratio can be related to the shape 
parameters as = M* and the minimum mass-to- 
light ratio is Tu{E) = '^mi{ehV — 100 for early type 
galaxies and T„(i) ~ 50 for late types in the (3 — 2 
model. 

According to our result the mass-to-light ratio of the 
brightest galaxies reaches much higher than 10'^. But 
this value is too high c ompared to those reported in lit- 
erature. For example, iMandelbaum et al] (|2006[ ) noted 
that a mass-to- light of early types reaches Mh/L^ = 674 
for their brightest samples (-22.5 < M < -22 in Tab. 4 
of their paper) of the SDSS lensing galaxies. The discrep- 
ancy may be explained by the fact that they are typically 
the brightest cluster galaxies. In our FSB halo identifi- 
cation method all mass not assigned to satellite galax- 
ies is assigned to the central subhalo after the bound- 
ness check, which makes the central subhalo as massive 
as the whole cluster and their mass-to-ligh t ratio very 
la rge (cf. Fig. 10 inlVale fc Ostrike^ |2006| and Fig. 4 
in [Tinker et al.l[2005h . Then, the group or cluster mass- 
to-light ratio is the one to be compared at the very high 
mass scales. 

4. ENVIRONMENTAL EFFECTS ON THE GALAXY 
DISTRIBUTION 




10" r 



108 ^ 



I I I 1 1 1 11 1 — I I I 1 1 1 11 1 — I I I 1 1 1 11 1 — r 



1 1 1 1 INI 


1 1 1 1 Mil 


1 1 1 1 Mil 


1— t 












= 2) Early 
= 2) Early fit 

{/c = 2) Early 

{/c = a) Late 



I I 



I I L 



10" 10'2 10'3 10" 

Mass (h-' Mq) 

Fig. 4. — The scaling relation between the halo mass and the 
galaxy luminosity (lower panel) or mass-to-light ratio (upper panel) 
in the = 2 and k = 2 model. In the lower panel, the solid 
line shows the derived luminosity of early type mock galaxies as a 
function of the subhalo mass and the short dashed line is the best- 
fitting curve for it. Also the early (long dashed) and late (dotted) 
luminosity distributions are shown for the k = 2 model. In the 
upper panel, the mass-to-light ratio of the mock galaxies is shown 
by the line of the same type as shown in the lower panel. But we 
do not plot the fitting result in upper panel. 

4.1. Definition of the Local Density 

For a comparative study of the environmental effects 
on the spatial and luminosity distributions of observed 
and simulated galaxies, a quantitative measure of the 
local environment is needed. To minimize the parame- 
terizations and maximize the spatial resolution we use 
the spline kernel to obtain the smooth galaxy number 
density field. It is given by 



Wiq) = 



i2-qf/{4nhl) 




for < g < 1, 
for 1 < g < 2, 
otherwise, 



where q = r/hs- The kernel is centrally weighted more 
than the Gaussian, and has only one parameter, hg. 
It is smooth to the second order and has a finite tail 
out to 2hs. Throughout this paper we set hg = d2o/2, 
where 6,20 is the distance to the twentieth nearest neigh- 
bor galaxy. Due to these features and to the adaptive 
nature of the resulting galaxy density map better 
represent the observed distribution of galaxies than the 
Gaussian, particularly in voids and clustered regions. 
To compare the observed and mock galaxy samples, it 
is necessary to make the number density of the den- 
sity tracers the same. We use galaxies brighter than 
= — 20 as the galaxy density tra cers. This selection 
differs slightly from that adopted by iPark et all (|2007al) 
who chose —21 < M < —20. We removed the bright cut 
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to better resolve the centers of clusters where galaxies 
brighter than are concentrated. 

4.2. Number Density Distributions 

The local density at galaxy positions is measured in our 
SDSS volume-limited samples. The density estimate is 
corrected for the boundary effects when the kernel sphere 
of radius d2o overlaps with the survey boundary. The 
late type galaxies having axis ratio b/a < 0.6 are re- 
moved from the sample to reduce the effects of internal 
absorption on luminosity. After this exclusion, we give 
the weight of 1/0.505 to late type galaxies to correct 
the galaxy number density for the missing inclined ones, 
where 0.505 is the fraction of late types with b/a > 0.6 
in a sample containing galaxies down to = —17.5 (the 
CM sample of Choi et al. 2007). 

The correction for boundary effect is unnecessary for 
mock galaxies since their distribution is periodic in all 
directions. The redshift distortions effects caused by the 
peculiar velocities are mimicked by shifting the mock 
galaxies along the x axis using x' = x + Vx/Hq, where 
is the a;-component of peculiar velocity and Hq is the 
bubble constant. 

Figure [5] shows the distribution of local density {open 
circles) at the location of galaxies in two volume-limited 
SDSS samples with absolute magnitude limits of —20.0 
(D5) and -18.5 (D2) for the = 2 model. Also plotted 
are the local density distributions of simulated galaxies 
selected by the same magnitude-limit criteria. It can be 
noted that the local density distributions of the SDSS 
galaxies and the corresponding mock galaxies match 
closely to each other at all densities except for at high 
densities. Also those distributions for the k = 2 model 
a re shown in Figure El with the same criteria. 

lOstriker et al.l (|2003f) employed a log normal function 
to fit the one-point distribution of dark matter den- 
sity and luminosity density in their hydro-simulation. A 
spherical top-hat filter of constant radius is adopted in 
their study. They found a good agreement between the 
simulated density distributions and the log normal distri- 
bution while the distribution of luminosity density shows 
a poor fit. We check how well the galaxy density distri- 
bution, instead of dark matter or luminosity, is described 
by the non-normalized log normal function 

"^^"dlogioA V2^a ' ^ ' 

where A = p/p and A is the amplitude of the distri- 
bution. The mean number density of galaxies is n = 
Alogj^o 6- The best-fitting results are shown with dashed 
lines in each panel of Figure [5] and the fitting param- 
eter values for various magnitude-limit (Alum) samples 
are listed in Table [H As can be noted, brighter galax- 
ies tend to be located at higher local densities. Also the 
local density distribution becomes broader for brighter 
subsamples if viewed in the linear scale of overdensity or 
has nearly a constant width {{a) ~ 1.29) in log-scale. 

Figure [7] and [8] shows the distribution of local density 
at locations of early (circles) and late (stars) type galax- 
ies in the (3 = 2 and k = 2 models, respectively. Open 
symbols are for the SDSS samples and the rest are for 
the simulated galaxies. The local density distributions 
of the simulated early and late type galaxies are very 
well matched with observations at low and intermediate 
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Fig. 5. — Distributions of local overdensities for the sample 
of mock {filled circles) and SDSS galaxies [open circles) in the 
magnitude ranges of M < —20 (upper) and M < —18.5 [lower 
panels) in the (3 = 2 model. The dashed line shows the log-normal 
fit to observations. 



TABLE 2 

The log normal function parameters of 
the local density distributions of the 
sdss galaxies 





sample 


A(/i-iMpc)-3 




cr 


-18.5 


D2 


4.53 X 10-2 


0.662 


1.27 


-20 


D5 


1.47 X 10-2 


0.949 


1.34 


-20.5 


D5 


7.44 X 10-3 


1.11 


1.30 


-21 


D5 


2.76 X 10-3 


1.23 


1.26 



densities {p/p < 10) even though there is only one input 
parameter in our model. Late type galaxies dominate 
these regions. From these figures, it appears that the 
galaxy morphology depends only on the halo mass and 
does not directly depend on environment. But the num- 
ber density analysis alone is not so sufficient to draw any 
definite conclusion of the environmental effect. We will 
further investigate this effect by comparing the luminos- 
ity functions of the SDSS with those of the mock galaxies 
in various environments in the next section. In the high 
density regions our model gives too few galaxies, which 
is probably due to the insufficient resolution of the simu- 
lation to maintain small subhalos within clusters and to 
the lack of gas physics. 
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Fig. 6. — Similar to Fig. [5]but for the k = 2 model. 

5. LUMINOSITY FUNCTION 

There have been several works reporting the environ- 
mental dependence of the galaxy luminosity function 
(Park ct al. 1994 and 2007a in observations and lMo et al.l 
[2004 in the HOD model). It was found that the char- 
acteristic galaxy luminosity (i*) is an increasing func- 
tion of the local density, and the faint-end slope (a) of 
the luminosity function is insensitive to the local density. 
However, Park et al. (2007a) used a spline kernel weight- 
ing to estima te local densities, and lMo eF al. (2004) (also 
lCooravll2005h use a spherical top-hat filter of a constant 
radius. For a quantitative comparison between observa- 
tions and models the same density estimation scheme is 
required. 

In the previous section, we studied the one-point dis- 
tribution of the local density at galaxy locations for each 
morphology sample of galaxies brighter than a certain 
absolute magnitude limit. In this section, we investigate 
the distribution of the absolute magnitude of early and 
late type galaxies located in different local density en- 
vironment. This gives us environment and morphology- 
specific luminosity function of galaxies. In our galaxy 
morphology assignment scheme we only use the observed 
total luminosity function and the parameter /3 (or k), 
the ratio of luminosity (or mass) of the early and late 
types at a fixed mass (or luminosity). It will be inter- 
esting to see how accurately our model reproduces the 
observed luminosity functions at different local densities 
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Fig. 7. — Local overdensity distributions for the two morpho- 
logical types in the (3 = 2 model. The distributions of local over- 
density of mock and SDSS galaxies (D5 and D2) are shown for 
two magnitude-limit samples M < —20 (upper) and A4 < —18.5 
(lower panels). Open cirlces and open stars mark the distribu- 
tions of overdensity of the early type and late type SDSS galaxies. 
And filled circles and skeletal stars are those of mock galaxies for 
the corresponding types. We tag each panel with the magnitude 
criteria of the sample in the upper-left corner. 



and for different morphological types. For this compar- 
ison we use the luminosity functions measured from the 
D3 sample of Park et al. (2007a) because this absolute 
magnitude-limited sample covers both bright and faint 
magnitudes well, relative to other volume-limited sam- 
ples. The fitting of the measured luminosity functions 
to the Schechter formula is carried out by using the MI- 
NUIT ^ packages which employ the maximum likelihood 
method. 

Figure [5] and [TO] show the comparison of the luminos- 
ity functions in four different local density regions in the 
P — 2 and k — 2 models. The resulting environment and 
morphology-specific luminosity function reproduces the 
observations surprisingly well. However, there are no- 
table disagreements for two cases. At very high densities 
{p/p > 10) the abundance of faint early type galaxies are 
significantly low in the simulation. This is again proba- 
bly due to the lack of small subhalos that were destroyed 
in high density environment. Another problem is seen 
at low densities {p/p < 1), where the early type galax- 
ies are too few in the bright end of luminosity function. 

^ http:/ /wwwasdoc. web. cern.ch/wwwasdoc/minuit/minmain. html 
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Fig. 8. — Same as Fig. [7]but for the k = 2 model. 

It seems that the morphology transformatin to bright 
early types in underdense regions through close interac- 
tions and mergers is more effi cient in the nature than in 
our model (jPark et al.ll2008l for observational evidence). 
In the high density regions, the k = 2 model describes 
the number distribution of faint late type galaxies better 
than the /? = 2 model while this k model shows a slight 
overestimation in the population of faint late types in 
the mean fields (0.4 < log(p/p) < 1) compared to the 
observation. 

Figure [Tl] and [12] compare the parameter of the 
Schechter function best fit to the SDSS data (open sym- 
bols) and simulated galaxy samples (filled symbols) as a 
function of local density in the /3 = 2 and k — 2 models. 
The dependence of both M.^* and a on local density is 
qualitatively well-reproduced by the simulation. How- 
ever, the characteristic magnitude Al* of the early types 
is significantly fainter at p/p < 2 in the simulation. This 
is due to paucity of bright early type galaxies in low den- 
sity regions, which is mentioned above. The parameter a 
of the simulation is quite different from the observation 
at p/p > 2 for early type galaxies. This is again due to 
the flat faint end slope of the luminosity function of the 
simulated early type galaxies. 

6. CENTRAL VELOCITY DISPERSION 

Figure [13] shows the central velocity dispersion of early 
type galaxies in the D3 sample. The gray symbols with 
connecting lines are the relations between tr^ and the lo- 
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Fig. 9. — Luminosity functions of mock {^curves) and SDSS 
(symbols) galaxies in the (3 = 2 model. Each line shows the lumi- 
nosity functions of the mock galaxies; solid lines for total, dotted 
lines for early type galaxies, and dashed lines for late type galax- 
ies. For comparison, we show the luminosity densities of the SDSS 
subsamples divided by the same density cuts. The sample selection 
criteria are based on the log scale of local density and are written 
in the lower-right corner of each panel. 
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FlG. 10. — Similar to Fig. [9]but for the ft = 2 model. 

cal density for early type galaxies in four subsamples with 
absolute magnitude limits of —19.0 > A4 > —19.3 (bot- 
tom curve). -19.3 > M> -20.1, -20.1 > > -20.5, 
and Ai < —20.5 (top curve with filled circles). The 
curves delineate the median value of cr„ in each local 
density bin, which monotonically increases as luminos- 
ity increases. We also plot the scaled one-dimensional 
velocity dispersion of the early type mock galaxies with 
black symbols. Because the observed velocity dispersion 
is obtained by sampling the inner part (< 1".5) of the 
galaxy, we have to apply a scaling factor to the mock 
galaxy velocity dispersions. Because it is irrelevant in 
this paper to derive the exact value of the scaling factor 
and apply it to the related analysis, we make a rough 
estimation of the relation between tr/os and T^, where 
aios is the line-of-sight aperture veloci ty dispersion and 
Vy is the virial velocity. According to iLokas fc MamonI 
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Fig. 11. — Dependence of the Schechter parameter values, 
{top) and a [bottom panel), on the local density for the SDSS D3 
{open symbols) and Mock {filled symbols) samples in the /3 = 2 
model. Circles mark the fitting values of full galaxy samples and 
boxes are for the early type samples. Open stars and filled hexagons 
are for the late types of the SDSS and mock samples, respectively. 
In this plot we add the error bars only to the SDSS fitting values 
and connect only the mock results with solid lines for clarity. 



(|2001f ). the resulting scaling ratio of the velocity disper- 
sions, /c(= cr;os/K), is 



0.5 < /c < 3 



(12) 



for the acceptable ranges of the velocity anisotropy and 
the concentration parameter in < (ra/Rv) < 1, where 
Va is the comoving aperture radius. Here we assume that 
Vy = ay, where ay is the three-dimensional velocity dis- 
persion and Vy = ^2GMy/Ry. In the range of /c written 
in equation (fT2|) . we simply set fc = 0.8 (M < —20.5), 
fc = 1 (-20.5 <M < -20.1), fc = 1.1 (-20.1 <M < 
-19.3), and fc = 1.2 (-19.3 < < -19.0) to match 
for the mean amplitudes of velocity dispersions for each 
mock and SDSS sample pair. At a fixed luminosity the 
velocity dispersion of the early type mock galaxies in- 
creases as the local density increases for bright galaxies 
and the slopes are decreasing for the faint sample. It 
is reassuring that this trend is exactly the kind of phe- 
nomenon found in the observation (gray symbols). While 
the switch of the slope in the faintest SDSS samples is not 
clearly observed in the mock sample (it shows a nearly 
flat slope on the local environments), the change of the 
slope is similar to each other. Both observation and our 
model show that the mass-to-light ratio of the early type 
galaxies is a function of environment, and that this de- 
pendence in turn is a function of luminosity. The mass- 
to-light ratio of the early types decreases as the local 
density increases for galaxies brighter than about Al*, 
but increases for those fainter than about A^* + 1. 
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Fig. 13. — Median distribution of central velocity dispersions of 
SDSS (gray symbols) galaxies and scaled one-dimensional velocity 
dispersions of Mock (black symbols) galaxies of early types in the 
four magnitude-limit samples divided by four magnitude cuts. The 
magnitude ranges of the subsamples are J\4 < —20.5, —20.5 < 
M < -20.1, -20.1 < M < -19.3, and -19.3 < M < -19 from 
top to bottom lines, respectively. 



7. SUMMARY & CONCLUSION 

We have proposed a model to assign galaxy luminos- 
ity and morphology to the dark subhalos directly iden- 
tified in cosmological N-body simulations. It is assumed 
that the galaxy luminosity is a monotonic function of its 
host halo mass. In the k model we assume that the halo 
masses of early and late type galaxies of equal luminos- 
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ity have a constant ratio. Another alternative model is 
the P model which assumes the constant luminosity ra- 
tio between equal-mass gala xies of different types . This 
model has been proposed bv lMarinoni fc HudsonI ([20M ) 
who found that the observed B-band luminosity func- 
tion of galaxies can be reproduced from the PS function 
assuming double power-law mass-to-light ratios and de- 
rived halo occupation numbe rs. It has been expanded by 
IVale fc Ostrikej (|2004l . l2006( l who adopted the satellite 
halo mass functions and directly link subhalos to the ob- 
served galaxies. In this paper, we have introduced the ra- 
tio of luminosity of the early and late morphological type 
galaxies at a given halo mass, and derived type-specific 
mass-to-light ratios as a function of subhalo mass. They 
are used to assign luminosity and morphology to subha- 
los. The mass-to-light ratio of the early type galaxies 
derived in this way has a minimum value of T„ ~ 100 at 
the scale of M„ ~ 3 x 10"/i"1Mq in the /3 = 2 model. 
The mass-to- light ratio starts to exponentionaly increase 
below Mu and increases in a power law above M„. 

We use the large-scale background galaxy number den- 
sity as an environmental parameter. The smooth galaxy 
density field is obtained by using the adaptive spline ker- 
nel which enabled us to resolve crowded regions well. The 
local density distribution of the SDSS galaxies is well de- 
scribed by the log normal function. The obtained log 
normal parameter values indicate that brighter galaxies 
tend to locate in dense regions. The local density distri- 
bution of the simulated galaxies is quite similar to that of 
the SDSS galaixes in voids and moderate-density regions 
for both morphological types. The underestimation of 
galaxy population in clusters is thought to be due to the 
evaporation of subhalos by the dynamical friction and 
tidal stripping. This can also explain the discrepance in 
the luminosity function of the early type galaxies in high 

density reg ions. 

Recently, iGott et all ()2006f ) measured the genus statis- 
tic from a large sample of the SDSS galaxies and 
compared it with those of mock galaxies created by 
the three distinctive methods such as the semi-analytic 
'alaxy form ation model applied to the Millennium run 
Springel et al.l 120051 ) . a hydrodynamic simulation, and 
the subhalo-galaxy correspondence model adopted in this 
paper. It was found that the observed topology of large 
scale structure was best reproduced by the subhalo- 
galaxy correspondence model even though other models 
were also consistent with the observation. However, the 
observed topology was marginally inconsistent with all 
simulations in the sense that it showed a strong meat- 
ball topology at the significance level of 2.5(7 at the scale 
studied. The prominence of the isolated high density re- 
gions in the observation seems to be due to the Sloan 
Great Wall which was a dominant structure in the sam- 



ple analyzed. 

We have found an impressive agreement between our 
simulated galaxies and the SDSS galaxies in the depen- 
dence of the central velocity dispersion on the local den- 
sity and luminosity. The early type galaxies tend to have 
higher or higher mass in high density regions at a given 
luminosity when they are brighter than about A^*. In 
other words, these bright galaxies tend to become rel- 
atively fainter in high density regions at a given halo 
mass. This interesting dependence of the mass-to-light 
ratio on environment was successfully reproduced by our 
subhalo-galaxy correspondence model of galaxy forma- 
tion. A more detailed study of this phenomenon will be 
presented in a forthcoming paper. 
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